Nucleation times in the 2D Ising model 
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A theoretical framework is presented for the estimation of nucleation times in systems with Brow- 
nian type dynamics. This framework is applied to a prototype system: the two-dimensional Ising 
model with spin-flip dynamics in an external magnetic field. Direct simulation results for the nucle- 
ation times, spanning more than four orders of magnitude, are compared with theoretical predictions. 
There are no free parameters; all system properties required by the theory are determined by simu- 
lations in equilibrium systems. Theoretical and simulation values are found to agree in most cases 
within 20%. 



I. INTRODUCTION 

Nucleation, the process of escape from a metastable 
state via the spontaneous growth of a stable nucleus in 
a metastable surrounding, has been studied often before, 
and excellent books and reviews exist [J. A prototype 
system to study nucleation phenomena is the well-known 
Ising model. Above the so-called critical temperature, 
in absence of an external magnetic field, up- and down- 
pointing spins are roughly equally abundant. Below the 
critical temperature, the system prefers to be in either 
of two states: one state with a positive magnetization 
in which most spins are pointing up, and one state with 
a negative magnetization. In the presence of an exter- 
nal field one of these states will become metastable, and 
will decay to the stable equilibrium state with a certain 
nucleation rate. 

In an earlier article [2|, we presented a theoretical 
framework to describe activated processes in systems 
with Brownian type dynamics. We then tested this 
framework in the study of magnetization reversal times 
in the two-dimensional Ising model. This first test-case 
was especially simple, since many properties of the tran- 
sition state, i.e. the state with highest free energy on the 
trajectory from one free energy minimum to the other, 
were well-known. For instance, it is a state in which half 
the system has a positive, and the other half a negative 
magnetization, and these two regions are separated by a 
pair of interfaces which span the shorter dimension of the 
simulation cell. 

Here, we apply this framework to nucleation phenom- 
ena. Starting in the metastable state in which most 
spins are anti-aligned with the external field, the system 
evolves via a transition state in which a droplet of spins 
aligns with the field, to the stable minimum in which 
most spins are aligned with the external field. This sys- 
tem is more complex than that of magnetization reversal 
in the absence of an external field: the size and shape 
of the critical droplet are not known a priori. Acharyya 
and Stauffer studied nucleation times in the two- threc- 
and four-dimensional Ising model with external field, and 
found agreement with prcdicitons of classical nucleation 
theory Here, we study the two-dimensional case in 
much more detail. 



The organization of our manuscript is as follows. In 
section^ we describe the model that we study in detail. 
Next in section lHll we outline the theoretical framework, 
based on our work on magnetization reversals, but ex- 
tended to nucleation phenomena. We then apply this 
framework to our prototypical model — nucleation in the 
Ising model with spin-flip dynamics. In section IIVI we 
compare the theoretical predictions with high-accuracy 
computational results. 



II. DETAILED DESCRIPTION OF THE MODEL 

We consider the Ising model on a square lattice with 
lateral dimension L, periodic (helical) boundary con- 
ditions, and an external magnetic field h which favors 
downward-pointing spins, with the Hamiltonian 



(1) 



in which = ±1 is the spin at site z, and J is the coupling 
constant. The summation runs over all pairs of nearest- 
neighbor sites; those of site i are j — i ±1 modulo N 
and j — i zt L modulo N, with N — LP' . The mag- 
netization is defined as M = ^ ^ Si ; it can take values 
M = —N, —N + 2, . . . , N; all through this manuscript, 
we restrict ourselves to systems in which L is even. As 
a consequence, M takes only even values, and summa- 
tions over a range of possible magnetizations only run 
over even numbers, with an increment of 2. 

The system evolves in time according to single-spin-flip 
dynamics with Metropolis acceptance probabilities 4] . If 
Ci is the configuration after i proposed spin flips, a trial 
configuration C^^i is generated by flipping a single spin 
at a random site. This trial configuration is then either 
accepted (Ci+i ~ C-^^) or rejected (Ci+i = Ci); the 
acceptance probability is given by 



Pn 



(2) 



in which /3 = l/{kBT) with Boltzmann constant ks and 
temperature T. The time scale is set such that in one unit 
of time, on average each spin is proposed to be flipped 
once. So in our system, in one unit of time we perform 
N Monte Carlo steps. 
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III. THEORETICAL FRAMEWORK 

To study the behavior of nucleation times at temper- 
atures below the critical one we may consider an ensem- 
ble of a large number of systems prepared in states with 
no large negative clusters present and study the rate at 
which one of the clusters in the system grows beyond the 
critical cluster size, after which the system is removed 
from the ensemble. The spin-flip dynamics described 
above may be represented by a master equation for the 
probability distribution P(S) of finding a system in the 
state S at time t. Due to the huge number of possible 
states this master equation cannot be solved analytically 
or even numerically for system sizes of practical interest. 
Therefore, as an approximation we assume that we may 
treat the clusters of negative spins in the system as be- 
ing independent. In this approximation, the probability 
that none of the Nc clusters has grown beyond the criti- 
cal size, is simply the iV*** power of the probability that a 
single cluster has not grown beyond the critical size. The 
dynamics of a single cluster is then modeled by a master 
equation for the equilibrium probability -P(C, t) that this 
cluster contains C spins at time t: 



dP{C,t) 
Jt 



Tc^c+iP{C + 1, t) + Tc^c-iP[C - 1, t) 



C+l,C 



Tc-i,c)P{C,t), 



(3) 



with Tec the transition rate from C to C . In section 
II V Bl we describe how to estimate the transition rates 

In order that the equilibrium distribution be a sta- 
tionary solution of the master equation we impose the 
condition of detailed balance 



r 



C+1,C 



exp [P{F{C)-F{C + l))], 



(4) 



where, up to a constant, (3F{C) — ~\nPeq{C), with 
Peq{C) the equilibrium probability of cluster size C. 

The long-time nucleation rate as predicted by the mas- 
ter equation Q follows as the largest eigenvalue (with 
a minus sign) of this equation, supplemented with an 
absorbing boundary at C = A. Here A is an integer 
larger than the critical cluster size Cx, chosen such that 
Peq{A) ^ Peq{Cx), and clustcrs with size A are almost 
certain to nucleate. The absorbing boundary condition 
is implemented by setting r^-i^^i equal to zero. 

The largest eigenvalue — of Tc,c' in Eq- ©; as well 
as the corresponding eigenvector Po(C), may be found by 
requiring that the net current away from cluster size C 
assumes the value i/Pq{C). Using conservation of proba- 
bility one easily checks that this may be expressed as 

jc+i,c = Tc+ucPoiC) -Tc,c+iPoiC + 1) 

= (5) 

c<C 

rA,A^lPo{A-l) 



(6) 



where jc+i.c is defined as the net current flowing from 
C to C + 1. This current may be approximated by 

_ Etiexp[-m^)] 
Jc+i,c - i^—c: . , c<Cx 

jc+i,c , c>Cx (7) 

because the sum on the right hand side of Eq. jSJ is dom- 
inated by the terms with small c- values, for which Pq{c) 
is approximately proportional to exp[— /3P(c)]. This 
may be checked in hindsight against the solution ob- 
tained. With this approximation the equation may be 
solved recursively for Po(c) in terms of Pq{A ~ 1) for 
c = A — 2, A — Z, . . . , with the result 



Po(c) _ Ta,a-i 



Po(A-l) 



A-l 



m+l,m " 



exp[/3(P(TO) - P(c))] 



Tm+l 



+ 1,171 



(8) 

Since the sum over m is dominated by values close to or 
larger than the critical cluster size, for which jm+i.m is 
approximately equal to we may replace Eq. |(SJ) by 



A,A-l 



^ exp[/3(P(m) - P(e))] 

m—c 



(9) 



Po(^ - 1) ^'^ ^ ^m.+l,7n 

Now substituting this into Eq. ^ we arrive at the result 



(10) 



where we have used the fact that the sum over c is dom- 
inated by small values of c to extend the sum over m to 
m ~1. The result in Eq. (|10|) is well-known. It is usually 
derived by considering a state with a stationary current 
in which mass is inserted at a constant rate on one side 
(e.g. at C = in our case) and taken out as soon as it 
reaches the absorbing boundary (see e.g. @, section IV 
E). In that case the replacement of jm+i,m by a constant 
is exact. 

One may give an even more accurate representation of 
the long time behavior of P(C, t) for arbitrary initial dis- 
tributions P(C, 0) that are concentrated near the origin, 
by writing it in the form 

P(C,t) = fcPo(C) exp[-vt] for t -> oo, (11) 

where fc is a constant which represents the overlap be- 
tween P(C, 0) and Po(C). This gives for the probability 
S{t) that the system has not yet nucleated at time t 

A-l 

S{t) = P{C, t) = exp[-u{t - td)] for t oo, (12) 

where td is called the delay time. The value of td may 
be obtained from the projection of the initial cluster size 
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distribution P(C, t) onto the most slowly decaying eigen- 
function of the master equation, (jSJ- This eigenfunction 
was given (with a different normalization) in Eq. 0: 



PoiA-l)rA,A-i 

A-l 



Po{n) 



^ . exp [piFjm) - F{n))] 



m+l,7n 



with im+i.m given in Eq. O. As a consequence of the 
condition of detailed balance, Eq. Q, the correspond- 
ing left eigenvector is obtained by multiplying ip^{n) by 
exp [/3F{n)], or 



A-l 



7U \ sr^ exp[/3i^(m)] 

V' (") = 2^ 3,n+l,m (14) 



Notice that for small values of n this is virtually indepen- 
dent of n, because only the largest m-values give impor- 
tant contributions to the sum. Using the proper normal- 
ization of the leading eigenfunction, and approximating 
P(C, 0) by 6c, 1, one now finds immediately that fol- 
lows from 



A-l 



^'(1)^ (exp[-/3F(n)]V?(n) 



exp [lytd] = 



71 = 1 



n=l 



[exp[-f3F{n)] (^'(n))' 



(15) 



Now, by subtracting unity on both sides, dividing by v, 
approximating sums over m from n to A — 1 by sums 
from 1 to ^ — 1, where appropriate, and using Eq. ^ 
to replace jm+i,m appearing in the left-hand sum of the 
nominator, one obtains 



W{n) 



miii[m,Ca:] 



A-l n-1 



J2 W{k) ^P\n) 



A-l 



^ (W^(n)^'(n)) 



n=l 



with W{n) = exp(-/3F(n)). 



(16) 



IV. SIMULATIONS AND RESULTS 

If one wants to apply the above theoretical framework 
to nucleation times in the Ising model, the two ingredi- 
ents required are: (i) the equilibrium probability Peq{C) 
to find a cluster of size C; and (ii) the transition rates 
rc',c from cluster size C to C". We obtain these two in- 
gredients via two different computational approaches. In 
all our simulations, we use a technique known as multi- 
spin coding , which enables us to reach long simulation 
times and thus good statistics. 



A. free energy landscape 

We measure the distribution of cluster sizes for vari- 
ous values of /3 and h, in a system with 64x64 spins. We 
are actually only interested in the distribution of clus- 
ters smaller than a certain size A, which was discussed 
in section IIIII Therefore we define a cut-off cluster size 
Cmax, which is typically chosen to be 300. We modify our 
algorithm in the following way: starting with a configu- 
ration Ci we perform a fixed number M of Monte Carlo 
steps, and measure the sizes of all the clusters of down- 
spins in the system. If there is a cluster with more than 
Cmax spins we reject the new configuration, and choose 
Ci-fi = Ci. Otherwise we accept the new configuration 
as our Ci+i. In both cases we add the distribution of 
cluster sizes in C^+i to the histogram. Finally we find 
the free energy F{C) from 



pFiC) 



In 



NiC) 



c<c„ 



NiC)' 



(17) 



where N{C) is the average number of clusters of size C 
in the histogram. 

If Cmax is chosen too large, the various clusters in the 
system infiuence each other. In particular, excluded vol- 
ume eff'ects around large clusters suppress bigger clusters 
in the metastable state more than smaller ones. As a 
consequence the distribution of cluster sizes depends on 
Cmax- To determine whether a certain value of Cmax 
is allowed, we check that for small clusters the free en- 
ergy curve coincides with the same curve, obtained with 
a lower value for Cmax- 

Figure n shows our measurements for the free energy 
according to Eq. (|17l) as a function of island size, for some 
combinations of temperature and external field. 

For large circular islands, classical nucleation theory 
can be used. The free energy is then approximated by 
the Becker-Doring expression 



F{C) w Fo + 2aV^ - 2hC, 



(18) 



where Fq is some constant, a the excess free energy of 
the interface per unit length, and h the strength of the 
external field. We fitted Fq and a to the data in figure ^ 
and added the corresponding curves as lines in the same 
figure. As expected, the measurements are well fitted by 
the curves, as long as the island size is not too small. 

The surface tensions obtained in this waj^are larger 
than those given by the Onsager expression 



a = 2J + /3"Mntanh/3J, 



(19) 



which is valid for long horizontal or vertical interfaces, 
and zero external field (see figure ^ . We will address 
this issue in a separate paper [9|. 
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FIG. 1: Free energy as a function of cluster size in the 64 x 64 
system at (3 J = 0.58 and h = 0.08 (o), (3 J = 0.53 and h = 0.08 
(+), and j3J = 0.56 and h — 0.06 (x) . In the simulations, C 
assumes integer values > 1. The lines represent the Becker- 
Doring expression Eq. II18II . with fitted values for a. 



1.2 




FIG. 2: Surface tension obtained by fitting the free energy 
curve, divided by the values given in Eq. 1191 . as a function 
of pj. Different symbols denote different strengths of the 
external field: h = 0.05(o), h = 0.06(0), h = 0.07(o), h = 
0.08(+), h = 0.09(«) and h = O.lO(x). 



B. interface diffusion coefficient 



The second ingredient in the theoretical framework in 
section Hill is the rate Tc+i,c of cluster growth. To esti- 
mate this rate we study the diffusion of a single interface 
in a system with anti-periodic boundary conditions in the 
absence of an external field, as described in 2,]. The loca- 
tion of the interface is obtained from the magnetization 



M. The diffusion coefficient D is defined as: 
■((M(i) -M(0))2\T 



D 



lim 

t — ^oo 



2t 



(20) 



This diffusion coefficient is, to a good approximation, 
found to be linear in the length of the interface B: 



D{B,L,l3J)^g{(3J)B 



(21) 



with the constant c very close to zero. The results for 
g{P,J) for various temperatures are plotted in figure|3| 
To arrive at an estimate for the rate Fc+i.c for clus- 
ter growth and shrinkage, we assume that the diffusion 
coefficient neither depends on the external field nor on 
the shape of the interface (straight or circular), but only 
on the length of the cluster boundary, which in our case 
we set to For large clusters this should be a very 

good approximation. From Eqs. I|1U|I and 1)16(1 we see that 
the contributions from small clusters to the expressions 
for v and td are almost insignificant, and inaccurate esti- 
mates of the cluster boundary length for those clusters is 
not very important. This then gives rise to a jump rate 



^C+i,C 



(22) 



The factor 4 arises because the jumps in magnetization go 
by units of 2. Notice that instead of identifying rc+i,c 
through D we could have taken Fcc-i-i just as well. Due 
to the free energy profile these quantities are not iden- 
tical, though never much different, but fortunately, for 
C-values close to the free energy maximum, which are 
weighted most strongly, their difference becomes very 
small. In addition there will be some compensating ef- 
fects, because for C > Cx one has Fc+i,c > ^c,c+i 
whereas for C < C'x this is just the other way around. 



C. nucleation rates 

To measure the long-time nucleation rates, we first 
bring the system into equilibrium. At time t = 0, we 
then reverse all the spins so that the system is near its 
metastable equilibrium. We measure the time the sys- 
tem needs to reach the magnetization corresponding to 
the stable equilibrium. We make a histogram of all the 
measured times. 

Figure 01 shows that at long times the decay function 
f{t) behaves as f{t) ~ exp(— i/^i). Here, we focus on 
the long-time behavior, and are specifically interested in 
the escape rate i>s- We obtain this quantity via a fitting 
procedure, in which we ignore the data up to a time tg, 
chosen such that f{t) shows exponential time behavior 
for t > to. Then we determine the time t' at which half 
of the remaining events have taken place. We then ob- 
tain the relaxation time r = from r = [t' — to)/ ln(2). 
Instead we could have made a linear fit of all the data 
points beyond to, but this makes no significant differ- 
ence. Table Ogives the measured relaxation times ti. In 
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FIG. 3; Monte Carlo measurements of the diffusion coefficient 
per interface length g{l3J), as a function of inverse tempera- 
ture pj . 
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5000 10000 15000 20000 

At 

FIG. 4: Histogram of the time /S.t elapsed before nucleation 
occurs, at inverse temperature /3J = 0.54 and external field 
strength h = 0.08, in a system with 64 x 64 sites. The straight 
line is a fit, given by n{At) ~ exp(Ai/1990). 



D. short-time behavior 

Besides the long-time exponential behavior of the nu- 
cleation probability, we also study the deviation from this 
behavior at short times. To do this, the master equation 
Q for the time evolution of one cluster is solved numer- 
ically, with the initial condition that the cluster consists 
of a single spin: 



(23) 



We then compute the cumulative probability distribution 
Pnuc{^,t) that the cluster has grown beyond a certain 
size Crnax during time t. The corresponding probability 
distribution for the nucleation of one of Nc statistically 
independent clusters at Cmax, at time t, is given by 



l-PnuciNc,t) = (1-P„„c(l,i)) 



(24) 



The quantity Pnuc{Nc, t) should be equal to the cumula- 
tive distribution of nucleation times, if we use as a def- 
inition of nucleation the first occurrence of a cluster of 
Cmax spins anywhere in the system. 

We have compared this result with the results of direct 
simulations, for j3 = 0.54 and h — 0.08. We did this 
for three different system sizes: 32 x 32, 64 x 64, and 
128 X 128. In all cases the starting configuration was a 
system with zero magnetization. Within a few time steps 
this develops into a quasi-equilibrium distribution as far 
as the distribution of small cluster sizes is concerned, so 
for the present system sizes the fact that we started from 
clusters of size zero rather than unity makes no difference. 
Figure [S] shows the results. Here, we used for N(, the 
mean number of clusters in the system as obtained from 
the same simulations that were used to obtain the free 
energy as a function of cluster size. 

The asymptotic slopes of the curves in the top panel 
of figure 13 correspond to the nucleation rates. The 
times at which straight-line fits to these curves cross 
1 — Pnuc = 1 correspond to the waiting times td- There is 
excellent agreement between the direct simulations and 
the parameter-free theoretical framework: The theoret- 
ical prediction obtained with Eq. ( I16|l is td = 234 MC 
time units, while the fits vary between td — 228 and 239 
MC time units. Also the behavior at short times is well 
described by the theoretical framework, as is evidenced 
in the bottom panel of figure [S] The agreement at late 
times was discussed above. 



V. DISCUSSION 



the same table, we also report the estimated relaxation 
times T2 = {{Nc) v)^^ with v as obtained in Eq. H10() and 
{Nc) = X]c'=i ^(C) the average number of clusters in 
the system. Finally we also give the ratio T2/T1. The 
table shows that, while the relaxation times span more 
than four decades in time, the estimated and measured 
relaxation times agree mostly within 20%. 



In this paper we showed that for two dimensional 
Ising models with spin-flip dynamics Bccker-Doring the- 
ory provides an excellent description of nucleation time 
distributions, provided a realistic description is used for 
the free energy of the growing droplets. We determined 
this free energy from cluster size distributions in equilib- 
rium Monte Carlo simulations and found it may be fit- 
ted well by the Becker-Doring expression, provided one 
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FIG. 5: 1 — Pnuc{i) obtained from the time evolution of one 
cluster (solid lines) according to the master equation JHJi com- 
bined with Eq. I|24|l . and by direct simulation of systems with 
32 X 32, 64 X 64, and 128 x 128 spins (the dotted lines) at 
PJ = 0.54 and h — 0.08. In both cases, the highest of the 
three curves corresponds to the smallest system, and the low- 
est one to the largest system. The lower plot shows the same 
data, but focuses on short times. 



From Eqs. and (^T)) it then follows that 



since the terms with m < mo in the first sum in l|10() 
basically do not contribute. As a consequence of this the 
asymptotic nucleation rate in the system is independent 
of the choice of mg. 

For the short-time nucleation rate the free energy of 
small clusters is also important, especially if one starts 
from a state of zero magnetization. However, if one starts 
from a quasi-equilibrium distribution representing the 
state after a sudden reversal of the magnetic field, the 
short time behavior will be dominated by clusters that 
were fairly large to start with and again Eq. (|17|l may be 
trusted. In 11] Van Beijeren gave explicit expressions for 
the short-time behavior, which may be used if the latter is 
well approximated by a diffusion equation in an external 
potential. In the present case these cannot be used, since 
on the relevant time scales the hopping process between 
neighboring cluster sizes is not well-approximated by a 
diffusion process. And since the hopping rates depend 
on cluster size no analytic expressions for the short-time 
behavior are available. But the numerical solution of the 
master equation Q gives a very good agreement with the 
results from our Monte Carlo simulation of the nucleation 
process, as was shown in Fig. |S1 

Besides the free energy as a function of cluster size 
our calculations require the transition rates Tc,c±i be- 
tween neighboring values of the cluster size. These we 
estimated by setting them proportional to the mean cir- 
cumference of a cluster, determining the proportionality 
constant from the simulated mobility of a straight inter- 
face in cylindrical geometry and imposing the detailed 
balance condition Q). 

In our estimations we have been using a number of 
assumptions, whose validity is not guaranteed under all 
conditions: 



uses surface tensions that are 10 to 20% higher than the 
surface tension of a bulk interface at zero field. Simi- 
lar conclusions were reached by Auer and Frenkelj^ in 
their studies of crystal nucleation in colloidal hard sphere 
systems. They also determined cluster free energies by 
monitoring the frequency of occurrence of clusters of a 
given size. Like we they could fit the resulting curves 
quite well by a Becker-Doring expression with an effec- 
tive surface tension that exceeds the bulk surface tension 
by in their case 20 to 40%. 

A point one may question is whether the assumption of 
independent noninteracting clusters holds even for clus- 
ters of small size. This may indeed be doubted, but for- 
tunately, at least for the asymptotic nucleation frequency 
(fNc) this is not really relevant. One may choose to de- 
fine as clusters only those clusters that have a size c > 
mo, with mo chosen such that clusters of this size already 
are very rare, but still exp(/3F(mo)) <C exp{(3F{Cx))- 



• Strong fields should modify the diffusion coefficient; 
this effect is neglected. The freedom to modify 
the field strength within the metastable region is 
limited though, and long nucleation times, as seen 
mostly in real experiments, require weak fields. 

• The diffusion coefficient is assumed to be deter- 
mined by the size of the cluster alone, and is calcu- 
lated on the assumption that this shape is strictly 
spherical. This requires that the temperature is not 
too low, because at very low temperatures the equi- 
librium shape of the cluster is squarish rather than 
circular [ij. On the other hand the temperature 
should not be too close to the critical temperature 
for shape fluctuations to be reasonably small. To 
some extent these fluctuations are taken into ac- 
count, since our calculation of the diffusion coeffi- 
cient is done for a fluctuating interface around a 
cylinder. But it is by no means certain that the 
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fluctuations of a circular interface arc in all aspects 
comparable to those of the interface around the 
cylinder. 

• It may happen that an island splits up or that 

two islands merge, corresponding in our theoreti- 
cal framework to non-zero transition rates Tc,c+i 
and Tc+i^c with i > 1. This effect also is partly 
accounted for through the numerical determination 
of the diffusion constant on a cylinder, but espe- 
cially for larger islands the difference in geometry 
may cause additional effects. However, these will 
only become important in systems that contain a 
sizable density of islands of spins aligned with the 
external field. Hence, also this approximation can 
be trusted least near the critical point. 

• No memory effects are accounted for explicitly. For 

spin-flip dynamics memory effects will chiefly be 
due to the influence of shape fluctuations on the 
transition rates. Since shape fluctuations on larger 
length scales will decay only slowly, these may be 
fairly long lasting effects. Again, our way of de- 
termining the difusion coefficient will take many 
of these effects into account implicitly, but mem- 
ory effects involving large shape fluctuations may 
be different for the present cluster geometry. For 
magnetization conserving dynamics (consisting e.g. 



of local spin exchanges) much stronger memory ef- 
fects exist due e.g. to the effect that a spin that is 
released from a large cluster has high probability of 
reattaching to it soon. 

Under conditions in which the effects above are neg- 

ligeable, our theoretical framework is able to estimate 
nucleation rates with an accuracy in the range of 20%. 
The small systematic ovcrcstimation by about 10% of 
the nucleation time by theory may have several causes. 
The radius of a cluster will be slightly larger than our 
estimate because especially a large cluster will typically 
contain a few holes in its interior. But in the investi- 
gated range of temperatures this should not amount to 
a correction of more than 3%. The assumption that the 
diffusion coefficient is independent of the orientation of 
the interface also may be not entirely correct. At low 
temperatures certainly a diagonal interface is much more 
mobile than a straight one, but at the present fairly el- 
evated temperatures one would not expect this effect to 
be large. Further there could be effects from the possi- 
bility of cluster splittings and mergings, though some of 
these certainly are accounted for through the numerical 
determination of the diffusion constant on a cylinder. 

We are currently extending our investigations to two- 
dimensional Ising systems with magnetization-conserving 
dynamics as well as to three-dimensional Ising systems. 
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TABLE I: Measured values (ti) and estimated values {T2) for 
the nucleation times in the 64 x 64 system. 
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